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INTRODUCTION 


Two projects are summarized in this report. The first project is entitled Stiffener Crippling ni- 
tiated by Delamination, and the second is entitled Pressure Pillowing of an Orthogonally Stiff- 
ened Cylindrical Shell. 


RESEARCH ACCOMPLISHED 


I. Stiffener Crippling Initiated by Delamination 

Background and objective 

Graphite-epoxy stiffeners tested in compression (Bonanni, Johnson, and Starnes, 1991) 
showed evidence of failure initiation in postbuckling by delamination at the free edges of the 
flanges. This mode of failure initiation occurred in /- and 7-section specimens with flange width to 
thickness ratios in the range from ten to twenty. The analyses conducted in the Bonanni study 
modeled the stiffeners as a branched shell using the STAGS computer code. Since shell theories 
arc based on plane stress and neglect through-the-thickness stress components, failure initiated by 
interlaminar stresses at the free edges of the flanges due to bending in postbuckling could not be 
studied with the STAGS model. 

The objective of this project is to develop a computational model of the stiffener specimens 
that includes the capability to predict the interlaminar stress response at the flange free edge in 
postbuckling. 

Presentation and Report 

A presentation of this work was given at a professional meeting. The citation for this presenta- 
tion is (speaker in boldface font) 

• Barlas, F. Aylin, and Johnson, Eric R., “Interlaminar Stress Concentration in Postbuckled 
Open-Section Stiffeners,” presented in session Nonlinear Structural Mechanics - 1, Joint 
SES/ASME/ASCE Meeting, University of Virginia, Charlottesville, Virginia, June 6-9, 
1993 

The following report was written and copies given to the technical monitor: 

• Barlas, F. Aylin, and Johnson, Eric R., “Variable Complexity Modeling of Postbuckled 
Stiffeners for Delamination Initiation,” Center for Composite Materials and Structures, 
Virginia Polytechnic Inst, and State Univ., Blacksburg, VA 24061, Rept. CCMS-93-1 1 and 
VPI-E-93-07, July 1993. 

This report was the basis of Ms. Barlas’s MS degree in Aerospace Engineering, which she 
successfully defended in June 1993. The abstract of this report follows. 

Delamination at the free edge is analyzed as a mode of failure for uniaxial compression of 
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postbuckled structural components. Analyses are performed for I-section stiffeners and for a 
dropped -ply laminate, all of which failed due to free edge delamination in earlier studies. These 
specimens were made ofAS4/3502 graphite-epoxy unidirectional tape. The nonlinear response of 
the specimens is modeled by a geometrically nonlinear finite element analysis. A variable-com- 
plexity modeling scheme is developed to economize computer resources for the nonlinear analy- 
sis. The three-dimensional stress field near the free edge is investigated by employing 20-node 
solid elements in that region. A dimensionless delamination index is calculated to investigate the 
severity of the interlaminar stress state. Numerical results correlated reasonably well with the 
experiments on the load-end shortening plots. The nodal lines and/or inflection points along the 
free edge in the postbucked configuration are found to be delamination critical sites. Both inter- 
laminar shear stress tangent to the free edge and tensile interlaminar normal stress are significant 
at these critical locations and are likely to initiate delamination. 

Interlaminar stress gradients 

An important finding in the latest work involves the mesh refinement necessary to capture the 
stress gradients in the vicinity of the free edge of the postbuckled flange. The gradient in the inter- 
laminar normal stress component in the direction of the normal to the edge is the fastest gradient, 
followed by the normal-direction gradient in the interlaminar shear stress component normal to 
the free edge, while the slowest gradients occur for the interlaminar shear stress component tan- 
gent to the free edge. Thus we found that the finite element mesh should be simultaneously 
refined in the direction normal to the free edge and through the thickness of flange in order to 
properly capture the stress gradient in the normal stress component. Mesh refinement did not sig- 
nificantly change the results for the interlaminar shear stress component tangent to the free edge. 
Early results obtained for specimens 12 and 110 were not converged to yield the results found in 
refined analysis in Sections 4.4.1 and 4.4.2 of the above cited report. 


II. Pressure Pillowing of an Orthogonally Stiffened Cylindrical Shell 

A conference paper was written on this project. The citation for it is 

• Johnson, Eric R., and Rastogi, Naveen, “Interacting Loads in an Orthogonally Stiffened 
Composite Cylindrical Shell,” Proceedings of Center for Adhesive and Sealant Science 
and Center for Composite Materials and Structures, Program Review, Donaldson Brown 
Hotel and Conference Center, Virginia Polytechnic Institute and State University, Blacks- 
burg Virginia, October 3-5, 1993, pp. II.M.l - II.M.15 

The text of this conference paper is given at the end of this report. 
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Interacting loads in an orthogonally stiffened composite cylindrical shell 


Eric R. Johnson 1 and Naveen Rastogi 2 

Aerospace and Ocean Engineering Department 
Virginia Polytechnic Institute and State University 
Blacksburg, VA 24061-0203 

Abstract 

The distributions of the interacting loads between orthogonal stiffeners and the shell wall are computed for a long 
circular cylindrical shell subjected to internal pressure. Identical stringers are equally spaced around the circumfer- 
ence, and identical frames, or rings, are equally spaced along the length. Both the stringers and frames are on the 
inside of the shell wall. Closed-end pressure vessel effects are included. A comparison is made between the linear 
e las tic response and geometrically nonlinear elastic response. The example problem has dimensions of the fuselage 
of a typical transport aircraft, and has a composite material shell wail and symmetric stiffeners. 

Introduction 

The cabin pressure in a transport aircraft causes about a 10 psi pressure differential across the skin. An unstiffened, or 
monocoque, fuselage would carry this internal pressure load as a shell in membrane response, like a pressure vessel. 
However, internal longitudinal and transverse stiffeners are necessary to carry maneuver loads, etc. The presence of 
these internal stiffeners prevents the fuselage skin from expanding as a membrane, and the skin bulges, or “pillows”, 
between the stiffeners under the action of the internal pressure. Where the skin is restrained against its expansion as a 
membrane along the stiffeners, a bending boundary layer is formed. Also, at the stiffener intersection a local concen- 
tration of the interacting loads between the stiffeners and the skin occurs. Hence, the objectives of this paper are 1) to 
analyze the concentration of the interacting loads at the stiffener intersection, and 2) to study the pillowing of the 
skin, for both a linear response and a geometrically nonlinear response of a fuselage subjected to internal pressure. 

Mathematical model 

An idealized mathematical model is assumed for the semi-monocoque fuselage to study the generic characteristics of 
the response in the vicinity of the stiffener intersection. The model is of a very long circular cylindrical shell inter- 
nally by identical stringers equally spaced around the circumference, and identical frames, or rings, equally 

spaced along the length. In general, the spacing of the stringers is not the same as the rings. The structure is periodic 
both longitudinally and circumferentially, and the loading is spatially uniform. Consequently, a structural repeating 
unit can be defined whose deformation determines the deformation of the entire structure. The repeating unit consists 
of a portion of the shell wall centered over portions of the stringer and ring as shown in Fig. 1. The radius of the mid- 
dle surface of the undeformed cylinder is denoted by R, and thickness of the cylindrical shell is denoted by t. Axial 
coordinate x and the circumferential angle 9 are lines of curvature on the middle surface, and the thickness coordinate 
is denoted by £, with -t/2 <,^<,1/2. The origin of the surface coordinates is centered over the stiffener intersection so 
that -l<,x<,l and - 0 ^ 9 £ 0, where 21 is the axial length, and 2RQ is the circumferential arc length, of the repeating 
unit 

The stiffeners are modeled as mathematically one-dimensional structural elements, such that the actions transmitted 
by the stiffeners to the inside surface of the shell wall are represented by distributed line loads as shown in Fig. 2. In 
this paper it is assumed that the stringer is symmetric about the x-C, plane through its centroidal axis, and the ring is 
symmetric about the &-C, plane through its circular reference axis. On the basis of the symmetry about the x- and 0- 
axes for the unit, only the interacting load components tangent and normal to the stiffeners are included in the analy- 
sis. The interacting load intensities per unit undeformed length along the contact lines are denoted by X^fx) for the 


1. Associate Professor 
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component tangent to the stringer, X^(x) for the component normal to the stringer, >^(9) for the component tangent 
to the ring, and by Xq(Q) for the component normal to the ring. The positive directions for these interacting loads act- 
ing on the shell are shown in Fig. 2. The purpose of the analysis is to determine these distributed line loads. 

An enclosed volume to contain the pressurized medium is modeled by assuming diaphragms to extend from the edges 
of the repeating unit to the axis of revolution of the cylinder. The diaphragms do not resist the deformation of the 
repeating unit, but act to transmit loads normal to the edges of the repeating unit due to the internal pressure. 

For both the linear and nonlinear response of the repeating unit to internal pressure, we use the Ritz method and the 
principle of virtual work augmented by Lagrange multipliers to enforce kinematic constraints between the structural 
components of the repeating unit The Lagrange multipliers represent the interacting line loads between the stiffeners 
and the shell. Displacements are separately assumed for the shell, stringer, and ring. 

Displacement and Lagrange multiplier approximations 

The periodic portions of the displacements are represented by truncated Fourier Series, and the nonperiodic portions 
of the displacements due to axial stretching are represented by simple terms in x. For the shell, the displacements of 
the middle surface are 


u(x, 9) = 

M N 

X I 

. „ % x 
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in which a = (mn) // and B = (nx) /© where m and n are non-negative integers. The displacements of the 

ftl ft 

centroidal axis of the stringer are 

M q x M 

«,(*) = X “rm sin + W r (x >= X W sm COS (0t m X > (4) 

«= 1 »=0 

Parameters q 0 in the axial displacement of the shell and q t in the axial displacement of the stringer represent elonga- 
tions of each respective element caused by either an axial mechanical load or due to closed-end pressure vessel 
effects. The displacements of the reference circle of the ring are 


N N 


M 0 ) = X V r« sin( Pn 9) 

w r ( 9) = X w ™ cos (M) 

(5) 

n = 1 

n - 0 


The distributions of the Lagrange multipliers, or interacting loads, are taken as 
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Note that in the stringer interacting normal load, the term A^q, which represents a uniformly distributed load, is omit- 
ted. Since the stringer as modeled by this analysis is not restrained from rigid body motion in the normal direction the 
resultant normal load on the stringer must vanish; i.e., 



( 8 ) 


\\ s {x)dx = 0, 


and this condition requires 10 vanish. 

Shell 

The generalized strain vector for the shell is 


Virtual work 


£ shell - [ £ xx’ e 90’ 7r0’ K xx’ K 00’ K J0^ 


(9) 


Assuming the strains are small and rotations are moderately small, the membrane strain-displacement relations are 1 

1 3v 
" R ( 30 


du 1.2 1 ,2 


. = -l— tv) ti^-V 


1 du dv 

Y *e ~ r%q + Tx + W »’ 


( 10 ) 


in which <t> x is the rotation of the normal around the 0-curve, <f>e is the rotation of the normal around the x-curve, and 0 
is the rotation around the normal. The rotations are given by 


1 <hv v . _ 1 f5v _ 1 du~ 

♦x _ “ Tx ~ Rde r * 2 Lax * ae 


cHv 


(ID 


The change in the normal curvature components in terms of the rotations are 


a^ 

K « dx 


1^0 


^0 1 * x 1 


K 96 “ RdQ dx + Rdd + R*' 


(12) 


The Donnell-Mushtari- Vlasov (DMV) approximation, or quasi-shallow shell theory, is obtained by neglecting the 
rotation about the normal in the strains of Eqs. (10) and (12), and the shaded term in rotation 4*9 °f Eq* (1 !)• 
assume that the contribution of the rotation about the normal to the strains is negligible for the shell, but do not 
assume the DMV approximation initially. 


Define the generalized stress vector in terms of the stress resultants and couples of Sander’s theory 1 " 2 by 


* shell 


- [**, ^00’ **0’ M xx< > 


such that the internal virtual work is 


8W ‘sHell ~ JJ & £ shetl shell dS ' 


(13) 


(H) 


where S denotes the area of the middle surface and dS = dx RdQ. For a laminated composite material shell wall, 
Hooke’s law is 


® shell ~ H£ shell 


H = 


A B 
B T D 


(15) 


-i 

in which the 3 x 3 submatrices A, B, and D are given by classical lamination theory . 
The statement of virtual work is 
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( 16 ) 


i 

5 Kheii = 8vv J If + J {*„ W »[« (x, 0) - (x, 0)] + \ s (x) 8w (x, 0) } dx 
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-e 


for every admissible variation in the shell displacements. For an enclosed volume, the work done by a constant hydro- 
static pressure is conservative. (Note that the pressure remains normal to the deformed surface.) In terms of shell dis- 
placements, the work functional for the constant hydrostatic pressure p is 


( 




( I + u.) {H- + (w 2 + v 2 ) ) + v y w + ^del 


U V w 

u x v x w x\ 

L u y v y 
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in which y = RB, and we denoted partial derivatives of the displacements by a subscript. The second term in Eq. (17) 
corresponds to the work of the pressure on diaphragms atx = ± / less the portion Q of this axial load that is carried by 
the stringer. The axial force Q is an additional Lagrange multiplier that accounts for axial load sharing between the 
stringer and shell. 


Stringer 

The virtual work statement for the stringer is 
l i 

5 { ( N xs de xs + M xs 5 k xs ) dx = - J { X xs W 5 1 U s (x) ~ e * W '* W 1 + W ^ W } ^ (18) 

-l -l 

+ Qq i 

in which Nn is the axial force in the stringer, M u is the bending moment, is the extensional normal strain of the 
r-gntrniHal line, is the change in curvature of the centroidal line, e s is the radial distance from the stringer centroid 
to the contact line along the shell inside surface, and the prime superscript denotes ordinary derivative with respect to 
x. The strain-displacement relations and Hooke’s law for the stringer are 


e xr = “'x + ^ w V 2 K x, = - w ’V 
N„ = (EA) s t xs M xs = (EI)k xs 


(19) 

( 20 ) 


Ring 

The virtual work statement for the ring is 

0 6 

J Whr + M Br 5 <Br^ R 0 " \ i\r^ V r + e rV + \r 5w r 

-e -e 

in which is the circumferential force, A/qt 1S *e bending moment, is the circumferential normal strain of the 
reference arc, is the change in curvature of the reference arc, and the radius of the contact circle between the ring 
and shell is R 0 + e r = R - 1/2. The rotation, strain-displacement relations, and Hooke’s law for the ring are 




( 21 ) 
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( 22 ) 


(23) 


♦, = ^ V r- W r> e 0r = ^<'V + M V> + 5^ **' = 

N 9 r= ( EA K £ er M *r= ( £/ >rV 
in which the over-dot denotes and an ordinary derivative with respect to 0. 

Displacement continuity 

Displacement continuity constraints along the contact lines between the stiffeners and the inside surface of the shell 


are 


8xs < x > = “ °) “ (jc ’ 0) " [“i W ~ e ^J\ = ° 

(*) = w (■*. 0) “ w s M = 0 

8&r (9) = v (0,0) -^ 9 (°,0) - [v r (0) + e,* r ( 0 ).] =0 

£<;,( 0) = w (0, 0) -w r (0) =0 


The variational form for these constraints are 
/ 


0 


(24) 

(25) 

(26) 

(27) 

(28) 
(29) 


j (S\s* X s + dx = 0 \ ( 5X 0r8e r + S \A ;r> + e r )dQ = 0. 

-I - 0 

The constraint that the elongation of the shell at 9 = 0 and the stringer are the same is 

8(2 [^i ~Q 0 ] = 0. 

Satisfying the variational form of the constraints, Eqs. (28) and (29), is not sufficient, in general, to satisfying the con- 
straints pointwise, Eqs. (24-27). Since the trigonometric Fourier Series functions are linearly independent, and Eq. 
(29) separately satisfies the polynomial terms in the axial displacements, the satisfaction of the variational constraints 
in this analysis does lead to satisfying the constraints pointwise. 

Discrete equations and their solution 

The discrete displacement vector for the shell is the (3MN +2M+2N +2) x 1 vector 

“shell = lit if. * 30) 

in which the subvectors are 

n r 

U Q = [<7o> W 00’ V 01’ W 01' *** V 0N* W 0N^ 

7 * 

“m = ^ u mO’ w mQ ,u ml’ v mV w ml , - - u mN ,v mN ,w mN^ m ~ ^ 

The (2M + /) x 1 discrete displacement vector for the stringer and the (2N + 1) x7 vector for the ring are 

“sir = [«1« “rl* w fl* ■’ u iM’ W smI T “r = v rV w rl v rN' w rN^' 

in which the term w s q for the stringer has been omitted since it does not deform the stringer and its conjugate resultant 
in Eq. (8) is equal to zero. This uniform normal displacement component is determined from the condition that the 
rigid body displacement of the stringer is the uniform portion of the normal displacement component of the shell at 0 

= 0; i.e.. 


(31) 


(32) 
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( 33 ) 


The discrete vectors of the Lagrange multipliers are 

Kr = [\c \sv X x*AT F ^ = l \r0’ X 9rV \rV — W V* 1 ^ (34) 

For the nonlinear problem an iterative solution procedure is required. At a fixed value of the pressure p, a sequence of 
displacements is defined by adding an increment to the previous member of the sequence to determine the next mem- 
ber in the sequence. For a good initial displacement estimate, the sequence converges to the displacement solution of 
the nonlinear problem. The initial estimate used here is the converged solution at the last pressure load step. The 
update procedure to determine the increment is based on the modified Newton-Raphson method. To obtain the equa- 
tions for the Newton-Raphson method, we first derived the virtual work functionals for the shell, stringer, and ring, 
based on replacing the displacement components by a generic member in their sequence plus an increment. Then the 
functionals are linearized in the increments. Second, the approximations in Eqs. (1) to (5) are substituted into the 
incremented virtual work functionals for the displacements, and similar forms are substituted for their increments. 
Integration over space is performed after substitution of Eqs. (6) and (7) for the Lagrange multipliers, and after substi- 
tution for the variational displacements and the variational Lagrange multipliers. (The variational displacements and 
the variational Lagrange multipliers use the same space of functions as their unvaried forms given in Eqs. (1*7).) This 
process results in a 3MN +6M + 6N +6 system of equations governing the increment in the displacements (indicated 
by a A preceding the displacement subvector symbol) and the Lagrange multipliers. These equations are 


*11 (*sheii)-P L ^shcli) 0 0 * 11 * 12*13 A U shell 

0 K 22 {u st r ) 0 B 2l 0 fl 23 Au slr 

0 0 K 33 (fl r ) 0 B 32 0 A u r 

b t u b t 2X 0 0 0 0 l str 

rL o bL ooo 3L 


Rshell ( u shell* P ) 
&str 

£r(«r> 


L <3 *23 0 0 0 L « J 

in which AT n , and A" 33 are the tangent stiffness matrices for the shell, stringer, and ring, respectively, that are 
functions of the displacements. The matrix L results from the nonlinear portion of the external work funcuonal for the 
hydrostatic pressure, Eq.(17). In the modified Newton-Raphson method these matrices K u , K 22 > K 33 , and L are only 
computed for the initial displacement in the sequence, and are not updated for each new member in the sequence. The 
submatrices By, i,j = I, 2, 3, in Eq. (34) are determined from the external virtual work terms involving the Lagrange 
multipliers, and the constraint Eqs. (28) and (29). These By submatrices are not functions of the displacements. The 
vector on the right-hand-side of Eq. (35) is the residual force vector, which vanishes at equilibrium for the repeating 
unit. The constraint equations correspond to the last three rows of the partitioned matrix in Eq. (35), and as shown in 
Eq. (35) these constraints are applied to the increments in the displacements. If the initial estimate of the displace- 
ment satisfies these constraint equations and the increments satisfy these same equations, then the final displacement 
in the sequence will satisfy these same equations. Equation (35) is solved first for the increments in the displacements 
in terms of the I ^grange multipliers, then this solution is substituted into the constraint equations to determine the 
Lagrange multipliers. Thus, the total solution is obtained. 

Numerical example 

Data for the example are 21 = 20 in., R = 1 17.5 in., 2RG = 5.8 in., f = 0.075 in., Rq = 113.72 in., e s = 1.10 in., e r — 3.78 
in., (EA) r = 0.592 x 10 7 lb., (EI) r = 0.269 x 10 8 lb-in 2 , (EA) S = 0.404 x 10 7 lb., (El), = 0.142 x 10 8 lb-in 2 , and with the 
shell wall stiffness matrices given by 


6 



lb-in . 



0.664 0.221 0 




"262 159 4.33" 

A = 

0.221 0.577 0 

_ 0 0 0.22 1_ 

*10 lb/in 

B = 0 

D = 

159 210 4.33 
4.33 4.33 159_ 


This data was originally used in an example by Wang and Hsu 4 , and the dimensional data is representative of a large 
transport fuselage All the results presented for this example are for an internal pressure p = 10 psi. The Fourier Series 
were truncated at M = N = 16, unless otherwise indicated. Since 20 = 2.83°, the shell in this example is shallow and 
the DMV shell theory should be adequate. We found that the numerical results using Sanders theory with the rotation 
about the normal neglected and the numerical results using DMV theory were essentially the same. 

Results and discussion 


Pillowing 

Circumferential distributions of the normal displacement w for the shell are shown in Fig. 3 for the linear analysis and 
in Fig. 4 for the nonlinear analysis. Axial distributions of ware shown in Fig. 5 for the Unear analysis and in Fig. 6 for 
the nonlinear analysis. For reference, the normal displacement for the unstiffened .shell, or membrane response, is w - 
0.2287 inches for the linear analysis, and w = 0.2290 inches for the nonlinear analysis. The presence of the stiffeners 
reduces the normal displacements from these membrane values as is shown in the figures. The w-distributions shown 
for the linear analysis compare very well with those presented by Wang and Hsu 4 . (Wang and Hsu’s results are lim- 
ited to linear analysis.). The pillowing effect is much more pronounced for the linear analysis (Figs. 3 and 5) than for 
the nonlinear analysis (Figs. 4 and 6). The largest normal displacement occurs midway between the stiffeners, and 
this value for the linear analysis is 0.1796 inches while it is 0.1541 inches when geometric nonlinearity is included. 
The minimum normal displacement occurs at the stiffener intersection (not shown in the figures), and its value is 
0.1392 inches in the linear analysis and 0.1490 inches in the nonlinear analysis. Normal displacements along the stiff- 
eners vary only slighdy from their values at the intersection for both analyses. Thus, including geometric nonlinearity 
in the analysis increases the minimum value of the normal displacement of the shell and decreases its maximum 
value, which is an indication that pillowing is reduced in the nonlinear response. 

The circ umf erential and axial normal strain distributions on the inner and outer surfaces of the shell also show the 
reduced pillowing effect in the nonlinear analysis. See Figs. 7-10. The circumferential bending strain (difference in 
eee between the inner and outer surfaces) is maximum at the stringer midway between the rings (Figs. 7 and 8), and 

the axial bending strain is maximum at the ring midway between the stringers (Figs. 9 and 10). These maximum 
bending strains are substantially reduced in the geometrically nonlinear response relative to the linear response. The 
values of the ^ and e^Q strains from the linear analysis (Figs. 7 and 9) compare very well to those presented by Wang 

and Hsu, except in one respect The exception is that the circumferential distribution of the axial strain at x = - 1 
(Fig. 9) does not exhibit a decrease in value as the stringer is approached. Wang and Hsu’s results, however, show 
decreasing to nearly zero as the stringer is approached along the circumference. We attempted several changes to the 
shell displacement approximations, and even programmed Wang and Hsu’s solution, but could not get any of these 
attempts to give a solution showing a decrease in the axial normal strain at the stringer. 

Interacting load distributions 

The distributions of the interacting loads between the ring and shell are shown in Figs. 1 1 and 12. The distribution of 
the circumferential component (Fig. 1 1) is antisymmetric about the origin, and has reduced magnitudes due to the 
geometrically nonlinear effect As shown in Fig. 12, the normal component of the interacting load is an extremum at 
die origin, exhibiting a severe gradient there. The negative value of at the origin indicates that the action of the 

ring is to pull the shell radially inward against the action of the pressure to expand the shell outward. The peak normal 
load intensity is reduced from - 1,674 lb/in in the linear analysis to -1,045 lh/in in the nonlinear analysis. 

The distributions of the interacting loads between the stringer and the shell are shown in Figs. 13 and 14. The distri- 
bution of the tangential component is antisymmetric about the origin and has reduced peak magnitudes due to the 
geometrically nonlinear effect as shown in Fig. 13. The normal component is maximum at the origin and has a 
steep gradient there as shown in Fig. 14. The maximum value of normal component is reduced from 484.7 lb/in in the 
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linear analysis to 320.3 lb/in the nonlinear analysis. Note that the resultant of this normal component vanishes as a 
result of Eq. (8). 

Stiffener actions 

The distribution of the force and moment resultants in the stiffeners are shown in Figs. 15 and 16. In the linear analy- 
sis the total axial load due to the 10 psi pressure is 3,407.5 lbs (thepK 2 © term multiplying q 0 in Eq. (17)), and the 
portion of this load carried by the stringer is 820.1 lbs, or about 24% of the total axial load, according to the value of 
Nj. at x = ± / shown in Fig. 15. For the nonlinear analysis at 10 psi, the axial load carried by the stringer is increased 

to about 861 .4 lbs. However, the axial force in the stringer increases slightly from its minimum value midway 
between the rings to a maximum at the ring location. The bending moment in the ring is reduced by the geometric 
nonlinearity at the expense of an increase in value of the circumferential force as shown in Fig. 16. The circum- 
ferential force in the ring is nearly uniform and increases from 7,242 lbs in the Unear analysis to 7,757 lbs in the non- 
Unear analysis. Thus, the force resultants in the stiffeners increase due to the inclusion of geometric nonUnearity in 
the analysis. 

Singularity at stiffener intersection 

The results for the interacting normal loads (Figs. 12 and 14) suggest that they exhibit singular behavior at the stiff- 
ener intersection. The argument for singular behavior is supported by plotting the magnitude of the interacting normal 
load at the stiffener intersection as M and N increase. The total normal load intensity acting on the cylinder at the stiff- 
ener intersection is denoted by where ^ = X^O) + X^(0). This value is -1,190 lb/in from the linear analysis and is 

reduced to -725 lb/in in the nonlinear analysis for M - N = 16. A normalized value of is plotted versus increasing 
M and N values, with M = N, in Fig. 17 for the linear analysis, and in Fig. 18 for the nonUnear analysis. The normal- 
ization factor, X^, is simply the value ofX^for the largest values of M and IV considered in each analysis. As shown 
in the figures, X^ is steadily increasing with an increasing number of terms in the truncated Fourier Series. Conse- 
quently, the series for X c does not exhibit, in the range of M and W considered, a convergent behavior. In spite of the 
fact that the normal load intensity at the stiffener intersection is exhibiting singular behavior, the total radial resultant 
load at the stiffener intersection converges rapidly with increasing values of M and N as shown in Figs. 17 and 1 8. The 
total radial resultant plotted in these figures is defined by 

e 

= J (X^ r cos8-X 0r sinfl) ( R 0 + e r )dQ , (36) 

-e 

since the resultant from the stringer vanishes by Eq. (8). From the linear analysis F ^ = -342.8 lbs, and from the non- 
linear analysis Fr = -375.5 lbs. Since the applied radial load due to internal pressure acting on the repeating unit is 
1,160 lbs (= 10 psi x 20 in x 5.8in), the ring resists about 29.6% of this applied load in the linear response, and this 
percentage is increased slightly to 32.4% in the nonlinear response. The remaining portion of the applied radial pres- 
sure load is carried by the shell. 


Conclusions 

The spatial distribution of the normal displacements of the cylindrical shell is more uniform, and the bending strains 
are reduced, in the geometrically nonlinear elastic analysis with respect to what is predicted by the linear elastic anal- 
ysis. That is, pillowing of the skin is reduced by inclusion of the geometric nonlinearity into the analysis. 

The series for the interacting normal load intensity at the stiffener intersection does not appear to converge in both the 
linear analysis and in the geometrically nonlinear analysis. The magnitude of the resultant normal load intensity at the 
stiffener intersection using a truncated series approximation is decreased by inclusion of the geometric nonlinearity 
into the analysis. However, the Fourier Series for the total radial resultant load carried by the stiffeners, which is 
resolved at the intersection, exhibits rapid convergence. This total radial resultant carried by the stiffeners is increased 
slightly by the inclusion of the geometric nonlinearity into the analysis. The axial force earned by the stringer due to 
the closed-end pressure vessel effect is increased in the nonlinear analysis with respect to its value in the linear analy- 
sis. Also, the circumferential force carried by the ring m increased in the nonlinear analysis with respect to its value in 
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the linear analysis. Thus, the stiffeners resist an increased portion of the internal pressure load, accompanied by a 
commensurate decrease in the load carried by the shell, when geometric nonlinearity is included into the analysis. 
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Fig.l Structural repeating unit of an orthogonally stiffened cylindrical shell subjected to 
internal pressure p. 



Fig.2 Distributed line load intensities (Ib/in) acting on the cylindrical shell due to the ring 
and stringer. 
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Fig.3 Circumferential distribution of the shell's 
normal displacement from the linear analysis. 



-1 - 0.8 - 0.6 - 0.4 - 0.2 0 

9/e 

Fig.4 Circumferential distribution of the shell's 
normal displacement from the nonlinear analysis. 



Fig.5 Axial distributions of the shell's normal 
displacement from the linear analysis. 


Fig.6 Axial distribution of the shell's normal 
displacement from the nonlinear analysis. 
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Flg.7 Circumferential normal strain on the inner and 
outer shell surfaces from the linear analysis. 
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Fig.8 Circumferential normal strain on the Inner and 
outer shell surfaces from the nonlinear analysis. 
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Fig.9 Axial normal strain on the inner and outer Fig.10 Axial normal strain on the inner and outer 

shell surfaces from the linear analysis. shell surfaces from the nonlinear analysis. 
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Fig. 13 Tangential component of the line load acting 
on the shell due to the stringer. 



X/l 

Fig. 14 Normal component of the line load acting 
on the shell due to the stringer. 









o 4 8 12 16 20 24 

M,N 

Fig. 17 Normal load Intensity A^and total normal load F^atthe stiffener 
intersection for increasing number of harmonics in the linear analysis. 
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Fig. 18 Normal load intensity A^and total normal load F^ at the stiffener 
intersection for increasing number of harmonics in nonlinear analysis. 





